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Abstract 

We have numerically studied the bosonic off-diagonal long range order, in- 
troduced by Read to describe the ordering in ideal quantum Hall states, for 
noninteracting electrons in random potentials confined to the lowest Landau 
level. We find that it also describes the ordering in disordered quantum Hall 
states: the proposed order parameter vanishes in the disordered (a xy = 0) 
phase and increases continuously from zero in the ordered (a xy = e 2 /h) phase. 
We study the scaling of the order parameter and find that it is consistent with 
that of the one-electron Green's function. 
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The quantum Hall effect (QHE) is a consequence of novel correlated states that arise 
in a two dimensional electron gas placed in a transverse magnetic field In ideal, i.e. 
translationally invariant, systems these states exist at isolated filling factors at which the 
system is incompressible. In systems with impurities these states broaden into phases - 
ranges of filling factor (v) which exhibit the same transport properties as the parent ideal 
states — thus giving rise to the characteristic plateaux structure of the QHE. The nature of 
the ordering in the ideal states was elucidated by Girvin and MacDonald (GM) and later 
by Zhang, Hansson and Kivelson (ZHK) || and by Read who showed that the Laugh- 
lin states could be viewed as condensates of composite bosons consisting of electrons that 
"carry" flux. The GM/ZHK formulation is distinct from that of Read and there is no proof 
that they are equivalent. In this paper we will be concerned, for reasons of computational 
convenience, solely with Read's formulation; we comment briefly on the GM case at the end. 

In this work we address the following questions: Is Read's bosonic off-diagonal long 
range order (ODLRO) a property of real, dirty quantum Hall systems, i.e. is it non-zero in 
the entire phase descended from an ideal state and does it vanish outside its boundaries? 
(That the ODLRO vanishes for a clean system when the Laughlin states are destabilized by 
varying the electron-electron interaction was shown already by Rezayi and Haldane || . ) We 
emphasize that at issue here is whether the ODLRO can serve as a sufficient characterization 
of the real systems that exhibit the QHE; note that incompressibility is lost when disorder 
is introduced into the system. Anticipating that the ODLRO does survive the introduction 
of disorder, we are led to the derivative questions of the critical behavior of the Read bosons 
near the transition between neighboring Hall plateaux, "seen" by them as a superfluid- 
insulator transition |J, and its relationship to the conventional measures of derealization. 

We will present likely answers to these questions based on numerical studies of non- 
interacting electrons subject to random potentials which are confined to the lowest Landau 
level (LLL) — a problem for which the localization properties of single-particle states have 
been investigated extensively J7|. This choice may seem perplexing on account of the his- 
torical association of the ODLRO with the fractional states, but we remind the reader that 
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the integer Hall state at v = 1, the filled LLL, is just the first state in the Laughlin sequence 
1, 1/3, 1/5... and does exhibit ODLRO; more generally, there is no distinction between the 
integer QHE and fractional QHE in this regard and we expect our qualitative results to 
hold quite generally. We note that even for non-interacting electrons the ODLRO is a prop- 
erty of the many-body state, and is not a one-electron quantity averaged over the occupied 
single-particle states. 

Most of the calculations are done on the sphere || using a density per flux quantum, 
Pi, of delta-function scatterers, and a modified version of Read's operator introduced by 
Rezayi and Haldane f5[. For a sphere containing = 25-1-1 flux quanta the single-particle 
states can be labeled by the eigenvalues, m, of L z which run from — S to S. If L z generates 
rotations about the point x, then Read's operator for the v = 1 state takes the form 

cf>Ux;S) = cl =s+1/2 F(x;S). (1) 

Here -F(x; S) is the flux insertion operator at x, which replaces a state of the system with 
flux quanta and occupations nj\^(m) by a state of the system with + 1 flux quanta 
and occupations riN^+iirn — 1/2) = n^,{m). Note that the electron creation operator acts 
on the states of the + 1 flux system. The intuitive content of this definition is more easily 
understood in the planar disk geometry, where the action of 4>r{0) is the insertion of one 
flux quantum at the origin causing each of the single-particle states to move outwards into 
its neighbor and the subsequent injection of an electron into the central Gaussian orbital. At 
v = 1 it is clear that this operation takes the N particle ground state to the N + 1 particle 
ground state, much as the operation of the field operator on the N particle ground state of 
a superfluid produces a state with an 0(1) overlap with the N + 1 particle ground state. 
For dirty systems, it is not immediately evident that (in obvious notation), (N + l\cj)\ i (x)\N) 
continues to be 0(1) for all v at which a xy = e 2 /h or that it vanishes when a xy = 0. 
There is, however, a suggestive connection of <fr R with the flux insertion in Laughlin's gauge 
argument [ IU| . The latter can be interpreted as the statement that flux insertion followed by 



transferring an electron between edges has no effect on the ground state in the v — 1 phase 



ljj , but it is important to note that in the gauge argument the electrons in the localized 
states are unaffected by adiabatic flux insertion while in the action of <j) R all the electrons 
are affected. We note that, in the gauge argument, the inertness of the localized electrons 
is essential in order to get a quantized a xy . Consequently, and in contrast, the expectation 
value of 4> R is not quantized. 

The principal numerical results concern the disorder-averaged absolute value of the an- 
tipodal correlator: 

G R {p,S) = [| (MO = 0)4(0 = *•))!]*., (2) 

where 6 = 0, 7r denote the north and south poles, respectively, and the expectation value () 
is evaluated in the iV e -electron ground state. (Note that the filling factor v = N e /N^.) We 
find: 

(1) In the thermodynamic limit (S — > oo), for v < 1/2, Gr{v, S) — > 0, while for v > 1/2, 
G R (v, S) — ► m R (i>) ^ 0, so that 4> R defines a legitimate order parameter for the dirty system. 
As seen from Fig. [I], m R vanishes continuously on approaching v = 1/2 from above. The 
critical point is at v — 1/2 because the random potentials employed were strictly particle- 
hole symmetric. We have also calculated the disorder- averaged absolute value of the order 
parameter matrix element, M R = [\(N e + 1|</4j(# = 0)|A^ e )|]di S , and checked that it behaves 
as expected. The antipodal correlator is, however, easier to calculate and vanishes more 
rapidly (exponentially with S rather than as 1/5") in the low-z/ phase, so it was investigated 
more thoroughly. For other practical reasons, we have focused on the sphere, rather than 
the disk or torus geometries. 

Previous work on destruction of the ODLRO in the ideal states p],£|,|5| identified the 
unbinding of zeroes of the wavefunction from the particles as the relevant mechanism. In 
our problem there is always one zero per particle due to fermi statistics; the ODLRO is 
destroyed instead by the derealization of extra zeroes (quasiholes) introduced by varying 
v. Heuristically, the variation of m R with v appears to reflect a geometrical property of the 
states, namely the extent of a "percolating" cluster of v — 1 liquid. The quantized Hall 



conductance, however, is sensitive only to the existence of the cluster. 

(2) For 1/2 < v < 1 one expects a state in which randomly localized quasiholes are 
interspersed with a v — 1 condensate. In the bosonic description, this state resembles a 
vortex glass where the phase of the order parameter at a given point in space fluctuates 



randomly between disorder configurations on account of fluctuations in the locations of the 
vortices (quasiholes). We have confirmed that not taking the absolute value in Eq. (||D causes 
the disorder average to vanish at all filling factors < v < 1. 

(3) The data in the critical region admit a finite-size scaling analysis consistent with 
the critical behavior found in previous studies of one-electron measures of localization. The 
appropriate scaling ansatz is 

G R ~S-^g[{v-\)S% (3) 

where x = l/(2z/^) and is the correlation length exponent. Note that the antipodal 
distance varies as the square root of 5*. We estimate > 2 and r\ ~ 1.6; consequently 
rn R {y) ~ [y — l/2) 2/3 where 2/3 = u^rj > 3.2. However, the analysis is complicated by a 
slow crossover, apparently reflecting a weakly irrelevant operator. The analysis is described 
below. 

(4) In order to further test the basic question of the order parameters for dirty quantum 
Hall systems, we calculated the antipodal correlations for the v = 1/3 order parameter 
(i.e., insertion of three flux quanta at a point followed creation of an electron there) in 
the noninteracting system. For finite S (and at p x = 4), the correlations are nonzero, with 
a maximum at a filling somewhat less than 1/3. The maximum value of the correlation 
function decays rapidly as the system size increases, perhaps as 1/S 3 . These results support 
the notion that there is a unique bosonic order parameter associated with each quantum 
Hall phase. 

We now turn to the details of the scaling analysis for Gr. We have data for Gr, as well 
as the equal-time one-electron Green's function, at p x = 2, 4, 8, 16 and for system sizes up 
to 2S = 320. At any given p ; , Eq. (§) describes the data well but the estimated (but not 



77) drifts substantially with the impurity concentration. 

Let us focus on the p\ = 4 data for the moment. The band-center data exhibit systematic 
deviations from scaling at the smaller system sizes, see Fig. 2; we estimate rj/2 = 0.78(2) 



based on the data for 25" > 120, but that value may be too low [ l3| . Likewise, a log- log plot 
of dGji/du at v = 1/2 (determined by fitting, as were the band-center values themselves) 
versus 2S exhibits curvature. However, the ratio of the slope to the value scales nicely, as 
shown in Fig. 2, leading to an estimate of x — 0.344(10), corresponding to = 1.45(4). 
The data collapse in the corresponding scaling plot, Fig. 3, is evident. This estimate of i/g 
above is significantly different from the value 2.34(4) determined by Huckestein and Kramer 
from the spatial decay of the Schrodinger Green's function. A similar analysis for the 
other impurity densities yields = 1.00(2), 1.64(5), 1.67(5) for p x = 2,8,16 respectively. 
These estimates suggest that p x is associated with a weakly irrelevant operator responsible 



for a slow crossover to the "quantum percolation" fixed point [14[, perhaps flowing from the 



classical percolation fixed point. In an attempt to investigate the issue of the correlation 
length exponent more closely we were led to calculations of the disorder-averaged absolute 
value of the antipodal Green's function, G e (^, S), for the same systems. Here we can only 
summarize the findings, which will be reported in detail elsewhere. At p x = 4 and for 
system sizes up to 2S = 241, we find that there is excellent data collapse (see Fig. 4) with 
z/g = 1.64(3). For the remaining impurity densities we estimate = 1.38(3), 2. 02(7), 2.16(8) 
for pi = 2,8,16 respectively. As the two sets of estimates for i/g, from Gr and G e , show the 
same systematics we do not view the difference between them as significant [[15] . 



A comment on the stability of our estimate for 77 is in order. This stability is somewhat 
surprising given the crossover effects in our estimates of v^. One possible explanation relies 
upon the numerically plausible identification of rj with the fractal dimension of the critical 
eigenstates which has been estimated to be approximately 1.6 |16|] and is not greatly different 
from the value 1.75 for classical percolation. If the latter is indeed the point of departure 
for the crossover the robustness of r] becomes plausible. 

To recapitulate, we have established that Read's operator does define a legitimate order 



parameter for the a xy = e 2 /h disordered quantum Hall system and by analogy, for other 
quantum Hall systems. At the transition between the ordered (a xy = e 2 /h) and disordered 
(&xy — 0) phases, we find that the ODLRO of the Read bosons exhibits a scaling behavior 
consistent with that of the one-electron properties. It is worth noting that the ordered phase 
is, in terms of the bosonic description, a superfluid state that resembles a vortex glass while 
the disordered phase resembles a bose glass; hence, there is a very suggestive analogy to the 



field tuned transition in dirty superconductors as discussed by Fisher |L7j in his treatment 



of the latter pfl . 

Finally, let us note some speculative implications for a couple of related issues. First, 
the consistent scaling of the ODLRO and Green's function implies, roughly speaking, that 
multiplication of the wavefunctions of the disordered Bose problem by the Jastrow factor 
T\ij( z i — Zj) is innocuous (at least as far as the correlation length exponent is concerned). 



This appears to give support to the claim of Jain, Kivelson and Trivedi [pl| that their 
"composite fermion" wavefunctions for the fractional QHE yield transitions in the same 
universality class as those for the integer QHE. Second, the GM/ZHK bosonization differs 
from Read's in that the former involves only the phase of the above Jastrow factor and that 
the resulting ODLRO is algebraic, even in the ideal states. Our cluster interpretation of 
the Read correlator results suggests that the GM correlator in the disordered phase will be 
characterized by the same v^. However, it is less clear whether the GM correlator exhibits 
the same algebraic decay thoughout the ordered phase, and we believe that this remains an 
interesting topic for future work. 
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FIGURES 

FIG. 1. The antipodal Read correlator, G R (is,S), plotted versus v for p\ = 4 at 2S = 20, 60, 
and 180 (from top to bottom). The data are averages over 3600 samples; the statistical errors are 
less than 1%. Inset: G R (v,S) plotted versus 2S for v = 0.9, 0.5, and 0.1, from top to bottom. 
Note that the v = 0.9 data saturates at large S, while at v = 0.1 it vanishes exponentially with 
S 1 / 2 . 

FIG. 2. Double logarithmic plots of G R (l/2, S) (x), dG R /dv(l/2, S) (O), and their ratio (+), 
versus 25. The line is a least-squares fit to the data for the ratio (2S = 40 through 320); the latter 
have been shifted vertically to accommodate them on the plot. 

FIG. 3. Scaling plot, of G R {u, S)/G R {l/2, S) versus S x (v - 1/2) at Pi = 4, for 2S = 20, 40, 

80, 100, 140, 180, 240, 320 and \u- 1/2| < 0.1. The continuous curve is the data at 2S = 320. The 
value of x used is 0.342, corresponding to = 1.46. 

FIG. 4. Scaling plot of the Green's function, G e (u, S), data at pi = 4 for all v and 2S = 41, 61, 

81, 101, 141, 181, 241. The ordinate is (ln[G e (l/2, S)/G e (y, S)]) 1 ' 2 and the abscissa is S x {v-\/2) 
with x = 0.31 corresponding to = 1.6. Note that G e (v, S) = G e (l — is, S). The different curves 
correspond to different S values. The data are plotted as curves rather than as points for clarity; 
the error bars are smaller than the width of the curves except for points near the origin. 
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